AUTOMATIC ENUMERATION OF REGULAR OBJECTS 



MARNI MISHNA 

Abstract. We describe a framework for systematic enumeration of families combinatorial structures which 
possess a certain regularity. More precisely, we describe how to obtain the differential equations satisfied 
by their generating series. These differential equations are then used to determine the initial terms in the 
counting sequence and for asymptotic analysis. The key tool is the scalar product for symmetric functions. 



Introduction 

Some classes of combinatorial objects naturally possess a substantial amount of symmetry and when formal 
sums of monomials encoding some parameter of interest are taken over the entire class, symmetric functions, 
or symmetric series appear. There has been some recent activity to determine how to extract enumerative 
series of sparse sub-families of these classes directly from the symmetric functions. The principle can be 
illustrated with one well studied example, the subset of labelled graphs in which the degree of each vertex is 
a fixed value, say k, known as the k-regular graphs. Here, we encode a graph by its degree sequence. When 
we consider the sum of this encoding over all graphs, they are encoded by the infinite product 

(1) G(X-L,X 2 , ...) =Y[0-+ x i x j)- 

i<j 

This is a well known symmetric series. Further, we remark that the coefficient of x\x\ ■ ■ ■ x\ in the formal 
power series development gives the number of labelled fc-regular graphs on n vertices. 

Such a coefficient extraction can be set up as a multidimensional Cauchy integral, as described by McKay 
for regular tournaments and Eulerian digraphs [20], and by McKay and Wormald for graphs with a fixed 
degree sequence |22j . However, in general this may not be a useful or practical formula. 

Indeed these techniques are well developed, and provide general formulae which we cannot currently obtain 
with the methods here, but they are difficult to make systematic, as they contain a saddle point analysis to 
make the asymptotic estimate which may be quite fine and specific to the problem. 

The primary goal here is to lucidly illustrate how techniques for computing the scalar product of symmetric 
functions in [7] can be a part of an essentially algorithmic process for asymptotic analysis. At the heart of 
the method is the fact that the scalar product of symmetric functions preserves a notion of D-finiteness [l4j , 
and, thanks to the algorithms in [7], this result is effective. 

We begin with a short recollection of symmetric series and D-finiteness, and a brief discussion on some 
places that D-finite symmetric series appear in combinatorics. We analyse graphs with fixed finite degree 
sets, and hypergraphs. Finally, in Section [4] we have the results of our semi-automated asymptotic analysis 
of these classes. 

1. Symmetric series and D-finiteness 

We provide a basic summary of symmetric functions in order to establish notation. The reader is directed 
to MacDonald's book [19] for full details. 

Denote by A = (Ai, . . . , Afc) a partition of the integer n. This means that n = Ai + • • • + A& and Ai > 
■ ' ' > Afe > 0, which we also denote A h n. Partitions serve as indices for the five principal symmetric 
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function families that we use: homogeneous (h\), power (p\), monomial (m\), elementary (e\), and Schur 
(s\). These are series in the infinite set of variables, Xi,X2, ■ ■ ■ over a field K of characteristic 0. When the 
indices are restricted to all partitions of the same positive integer n, any of the five families forms a basis for 
the vector space of symmetric polynomials of degree n in x\, X2, ■ ■ ■ On the other hand, the family of p^s 
indexed by the integers i € N generates the algebra A of symmetric functions over K: A = K[p\,p2, ■ ■ ■]■ 
Furthermore, the p^s are algebraically independent over Z. 

Generating series of symmetric functions live in the larger ring of symmetric series, if [t][[pi,_P2j • ••]]• 
There, we have the generating series of homogeneous and elementary functions: 

H(t) = Y J h n t n = exp (X>i) ' E ® = E e « tn = ex P (E(- 1 )>4 J ■ 

n \ i / n \ i / 

We often refer to H = H(l) and E = E(l). 

Alternatively, the power notation A = l" 1 • • • k nk for partitions indicates that i occurs rii times in A, 
for i = 1, 2, . . . , k. The normalization constant 

z x := r i n 1 l---k n "n k \ 

plays the role of the square of a norm of p\ in the following important formula: 

(2) (p\,P») = 5\^z x , 

where is 1 if A = /j, and otherwise. 

The scalar product is a basic tool for coefficient extraction. Indeed, if we write F{x\,X2, ■ ■ ■) in the 
form ^ A fxmx, then the coefficient of a;^ 1 • • • x k k in F is fx = (F, h\). Moreover, when A = 1™, the identity 
h\n — pin yields a simple way to compute this coefficient when F is written in the basis of the p's. When 
viewed at the level of generating series, this fact gives the following theorem: 

Theorem 1.1 (Gcsscl[14 ; Goulden & Jackson[15j). Let 8 be the K-algebra homomorphism from the algebra 
of symmetric functions over K to the algebra K[[t\] of formal power series in t defined by 9{p\) = t, 9{p n ) = 
for n > 1 . Then if F is a symmetric function, 

n=0 

where a n is the coefficient of F. 

To end our brief recollections of symmetric functions recall that plethysm is a way to compose symmetric 
functions. An inner law of A, denoted u[v] for u, v in A, it satisfies the following rules [29], with u, v, w G A 
and a, in K 

(au + 0v)[w] — au[w] + /3v[w], (uv)[w] = u[w]v[w], 

and if w = ^2 X c x p x then p n [w] — J2x c \P(nX 1 )P(nX 2 ) ■ ■ ■■ For example, consider that w[p n ] = p n [w], and in 
particular that p n [p m ] — Pnm- I n a mnemonic way: 

w\p n ] =w(pin,P2n,- ■■,Pkn,---) whenever w = w(pi,p2, ■ ■ ■ ,Pk, ■ ■ •)■ 

1.1. D-finite multivariate series. Recall that a series F £ K^X\, . . . , x n ]] is D-finite in xi, . . . ,x n when 
the set of all partial derivatives and their iterates, <9 nH Hn F/dx 1 ^ ■ ■ ■ dx l ™ , spans a finite-dimensional vector 
space over the field K {x\, . . . , x n ). A D-finite description of a series F is a set of differential equations which 
establishes this property. A typical example of such a set is a system of n differential equations of the form 

df d k f 
1i{x)f(x) + q 2 {x) — (x) H h q k (x)—j;{x) = 0, 

where i ranges over 1, . . . , n, each qj is in K(x\, . . . , x n ) for 1 < j < k, and k and qj depend on i. 

Such a system is a typical example of a D-finite description of a functions, and often this will be the 
preferred form for manipulating /. In truth we can accept any basis which generates the vector space of 
partial derivatives, but in the applications below, this form is particularly easy to obtain. 



1.2. D-finite symmetric series. The following definition of D-finiteness of scries in an infinite number 
of variables is given by Gessel [14] , who had symmetric functions in mind. A series F £ K[[x%, X2, ■ ■ •]] is 
D-finite in the Xi if the specialization to of all but a finite (arbritrary) choice, S, of the variable set results 
in a D-finite function (in the finite sense). In this case, many of the properties of the finite multivariate case 
hold true. One exception is closure under algebraic substitution, which requires additional hypotheses. 

The definition is then tailored to symmetric series by considering the algebra of symmetric series as 
generated over K by the set {pi,p2, ■ • ■ }: a symmetric series is called D-finite when it is D-finite in the PiQ. 

Example. Both H{t) and E(t) are D-finite symmetric functions, as for any specialization 
of all but a finite number of the p^s to results in an exponential of a polynomial. Similarly, 
exp(/ifci) is D-finite because hf. = ^2\\- k P\ is a polynomial in the piS. 

The closure under Hadamard product of D-finite series |18j yields the consequence: 

Theorem 1.2 (Gessel). Let f and g be elements of K \t\, . . . , ife][[pi,p2i • • •]]; D-finite in the pi's andtj's, 
and suppose that g involves only finitely many of the pi 's. Then (/, g) is D-finite in the tj 's provided it is 
well-defined as a power series. 

1.3. Effective calculation and algorithms. In our initial study [7] we gave an algorithm which, given 
a D-finite descriptions of two functions satisfying the hypothesis of Theorem 11.21 determines a D-finite 
description of the series of the scalar product. Henceforth, we shall refer to this algorithm as scalar_de. 
As we noted in [7], a second algorithm, hammond, based on the work of Goulden, Jackson and Reilly [15] 
applies in the case when g = exp(/i„i), which we shall sec is precisely how one can extract the exponential 
generating series of sub-classes with "regularity" . They are implemented in Maple, are are available for 
public distribution at the website http://www.math.sfu.ca/~mmishna. Maple worksheets illustrating the 
calculations presented are also available at that same site. 



Species theory (in the sense of [3j[T7]) is formalism for defining and manipulating combinatorial structures 
which relates classes to encoding series. An important connection to our work here is that the series for 
structures we consider are D-finite symmetric series, and many of the natural combinatorial actions preserve 
D-finiteness on the level of these series. 

The reader unfamiliar with species is heartily encouraged to consult [3]. A species associates to every 
set a family of structures in a way such that two sets of the same cardinality yield the same family, upto 
isomorphism. For example, the species of sets E on the underlying set U is simply E[[/] = U. The species of 
lists \-[U) — {(xi, X2, ■ ■ ■ , x n ) : Xi € U, n = 0, 1, 2 . . . }, is the set of finite ordered collections of elements. The 
atomic species, X[U] is U if U contains a single element, and is empty otherwise. 

The theory of species develops a rigorous formalism which allows a sort of calculus of combinatorial 
families. For example we construct lists of length 4 from our atomic species via multiplication: \-n[U] = 



The key feature that we use are that for every combinatorial family (species) F that one can define, there 
is an associated cycle index series Zp and an asymmetric cycle index series Tf both of which are symmetric 
series. Recall for any species F its cycle index series Zf is the series in C[pi,p2, ■ ■ ■] given by 



where the value of Fix F[A] is the number of structures of F which remain fixed under some labelling permu- 
tation of typeQ A, and mu gives the number of parts of A equal to k. 

The definition of the asymmetry index series of a species F, denoted Tf, as introduced by Labelle [3] is 
related, but more subtle. The series T behaves analytically in much the same way as the cycle index series, 
notably, substitution (in almost all cases) is reflected by plethysm, etc. Essentially, this series counts the 
objects with no internal symmetry. Table [1] contains some small examples of both series. 



This is interestingly enough not equivalent to D-finiteness with respect to either the h or e basis. 
2 A permutation of type (l m i , 2 1 ™ 2 , . . .) has mi fixed points, ni2 cycles of length 2, etc. 
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2. D-FINITE SYMMETRIC SERIES APPEAR NATURALLY IN COMBINATORICS 
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Table 1 . Index series of small species and their corresponding symmetric functions 



In a fashion similar to the cycle index series, r F arises through the enumeration of colourings of asymmetric 
F-structures. 

A notable example is the species of sets, E. Recall for any finite set U we have that E[U] = U. The two 
series above turn out to be Z E = expQ^ n p„/n) = J2n h n and T E = expQ^ n (-l) n p n /n) = J2 n e «- 

The primary advantage of this approach, as is true with any generating series approach, is that nat- 
ural combinatorial operations (set, cartesian product, substitution) coincide with straighforward analytic 
operations (sum, product, plethystic substitution 

The exponential generating series of a species F is the sum F(t) = J^ ra |F[n]|^j, where |F[n]| is the number 

of structures of type F on a set of size n. The ordinary generating function, F(i), is the sum F(t) = 
J2 n Orb(F[n])f", where Orb(F[n]) is the number structures of F on a set of size n distinct up to relabelling. 
Also recall the notation [x n ]f(x) refers to the coefficient of [x n ] in the expansion of f(x). This definition 
extends likewise to monomials. 

The next result is essentially a collection of known results and basic facts of D-finite series. 

Theorem 2.1. Suppose F is a species such that Zf is a D-finite symmetric series and write p n = X1+X2 +. . .. 
Then all of the following series are D-finite with respect to t: 

(1) The exponential generating function F(i); 

(2) The ordinary generating function F(t), if the additional condition that Zf(pi,p2, ■ ■ ■) is D-finite with 
respect to the X4 variables is also true; 

(3) The series J2 n {i x i ' ' ' x u} Z f) for fixed k; 

(4) The series J2 n Efc G s» ( x i lx 2 ■■■ x n n z ^ t n /n\, for any finite set ScN. 
Proof. The first two parts are proved using two basic results about cycle index series: 

F(i) = Z F (t,0,0, ...) and F(t) = Z F (t, i 2 , i 3 , . . .) 

The first specialization is well-known to preserve D-finiteness [35] for any n. The additional condition on the 
second item is sufficient to prove the D-finiteness since the stated substitution is the same as £1 1— > t, and 
Xi 1 ► otherwise. 

The third item of the proposition is proved by the expression 

in 

^([^••• 2 ; I fc jZ F )- = (Z F ,exp(^)), 

n 

which is D-finite by Theorem 1 1.21 

The final item of the proposition is true because the series is equal to 



Z F ,expfty^ hj) 



'This is slightly less true with the asymmetry index series, but true enough for our purposes. 
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Figure 1. The graph associated to the coloured set partition {{I3, 42}, {22, 64}, {34, 73}, {52, 81}} 

which is also D-finite by Theorem 11.21 □ 

We have one large class of species for which the cycle index series is D-finite. All of our examples come 
from this class. 

Theorem 2.2. Let E be the species of sets and let P be a polynomial species with finite support. Then F = 
E o P describes a species for which Zf is a D-finite symmetric series and provided that P(0) = 0, Tf is also 
a D-finite symmetric series. 

Proof. If P is a polynomial species, then its cycle index series is a polynomial in the Pi's, say P(pi, ■ ■ ■ ,p n )- 
Composition of species is reflected in the cycle index series by plethysm, thus 

Z F = exp(^2p k /k)[P(pi, . . . ,p n )] = exp {^P(p k ,p 2k , ■ ■ ■ ,Pnk)/k) ■ 

k 

For any specialization of all but a finite number of pi to zero, this gives an exponential of a polynomial, 
which is clearly D-finite. Thus, Zf is D-finite. We can similarly show that Tf is also D-finite under the 
stated conditions, since the composition also results in a plethystic composition. □ 

Our concluding remarks in Section 15.21 address the more general question of combinatorial criteria on a 
species F that ensure that Zf or Tf are D-finite. 

3. Using species to describe regular graph-like structures 

Ultimately our goal is to generalize the well-studied case of fc-regular graphs to other structures whose 
cycle index series are D-finite. To do so, we express the graph encoding by degree sequence as symmetric 
series, and describe how to find such a representation in general using species theory. 

In Eq. ([!} we define G{x\, X2, . . .) as the encoding over all graphs of their degree sequence and we express 
this as an infinite product. It turns out that this series is equivalent to E\e^\ , which is equivalent to Fe e 2 • The 
equivalent series for multigraphs (with loops) is equal to H[h2\ — Zf Q f 2 , and thus suspecting an explanation 
via species, we investigate this connection. Specifically, how do we construct symmetric function equations 
to describe the generating functions of different families of objects, such as hypergraphs. 

We begin with the remark that F = E o E2 is not the species of graphs. It is the species of partitions into 
2-sets. For example, {{1, 4}, {2, 6}, {3, 7}, {5, 8}} is an element of F, and we should not think that this is the 
graph on 8 vertices, with four edges, rather it just gives the basic structure, i.e. four edges. 

We express the Polya cycle index in the power series symmetric function. As a series in the symmetric 
Xi indeterminates, it is an inventory of distinct (non-isomorphic) colourings of the elements of the species. 
For example, the non-isomorphic colourings (by positive integers, say) of the set {a, b} is the set of maps 
{(a, b) 1 !■ € N 2 : i < j}, and the inventory of all such colourings is J2i<j x i x j = ^2- 

A colouring of an element in E o E2 gives rise to a graph (See Figure [T]) and two colourings are isomorphic 
if one is a graph relabelling of the other. The monomial encoding a colouring indicates how many time each 
colour was used, that is, in how many edges the colour appears, that is, the degree of the vertex represented 
by the colour. 

We restate this correspondence. The species E o E2 indicates the structure- sets of pairs. The cycle 
index series ^eoE 2 encodes non-isomorphic colourings of elements, which are in turn equivalent to labelled 
multi-graphs. The sets of pairs indicate edges, and the colours indicate vertices. 




Species which can be written as polynomials of molecular species. For example, every species in Table [T] is polynomial. 
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Figure 2. A structure composed of a set of smaller structures (cycles and 3-sets). 



For many applications, like regular graphs, we would like to count colourings without repetition. In this 
case, we do not allow repetition of a colour in a given object, hence to encode a fc-set, each colour appears 
exactly once, and this is precisely the notion of asymmetry in the asymmetry cycle index, and thus we use 
r instead of Z . Remark, and thus 

^E k = <Sfc and thus, = E. 

Taking the same species E o E2 as above, and using the asymmetry index series with a similar argument, we 
get that TEoEa = E\e-i\ encodes simple graphs without loops on the set of colours precisely as is determined 
by Eq. ([1]): E[e2\ — Yii<jO- + x i x j)- This gives us a way to have direct access to monomial encodings of 
combinatorial objects, as symmetric functions expressed in common bases, like the power sum basis. These 
two series are compatible and, one can show that graphs with loops are encoded by E[h,2], and graphs with 
multiple edges, but no loops are given by H[e 2 ]. 

More generally, we can consider any structure which is built as a set of objects from a finite set of classes. 
Figure [2] shows a more general object built as a set of cycles and sets. In this framework it is encoded by the 
monomial g, and thus we see that regularity in this situation refers to the number of times each 

label appears in one of the smaller substructures. 

Using this framework we can examine other species of structures built up from smaller objects. These 
species are such that both Zf and Tf give rise to interesting combinatorial objects. 

We can produce enumerative results for objects all of the same flavour: labelled sets of objects in which 
there is a certain regularity. We begin with a natural generalization of fc-regular graphs, and then we consider 
other types of objects such as hyper-graphs, and directed graphs. 

3.1. S- regular graphs. A graph is S -regular if the set of vertex degrees in the graph is a subset of S. For 
example, a graph is {i, j}-regular if every vertex is of degree i or j. 

It does not seem that the asymptotic enumeration of these objects has been directly considered before. 
It is, in some sense, a variation of the asymptotic number of labelled graphs with a given degree sequence, 
which has been considered by Bender and Canfield [2] and McKay and Wormald [22] , and may very well be 
computable from this. 

Thus, the scalar product which represents the generating series for the number of {i, j}-regular graphs 
is given by incorporating this factor, which ultimately greatly simplifies the calculation. The exponential 
generating series for the number of {i, j}-regular graphs is given by 

G itj (t) = (£[e 2 ],exp {t(hi + hj))). 

This is clearly D-finite, and computable using scalar_de (although not hammond). Furthermore, by a 
similar computation, we have the following result. 

Theorem 3.1. The number of S -regular graphs is D-finite for any finite ScN, and its exponential gener- 
ating series is given by the scalar product of symmetric functions, 



i, j Initial terms in the counting sequence 

1.2 1,0,1,4,18,112,820,6912,66178,708256,8372754,108306280,1521077404 

1.3 1, 0, 1, 0, 8, 0, 730, 0, 188790, 0, 102737670, 0, 102172297920,0 

1.4 1,0,1,0,3,6,30,1011,38920,1920348,116400186,8580463110,757574641296 

2.3 1,0,0,1,10,112,1760,35150,848932,24243520,805036704,30649435140 

2.4 1,0,0,1,3,38,730,20670,781578,37885204,2289786624,168879532980 

3,4 1,0,0,0,1,26,820,35150,1944530,133948836,11234051976,1127512146540 

Table 2. Counting sequences for {^jj-regular graphs for small values of i and j. 



Table [2] offers the initial counting sequence for some small values of fc and j . Each one corresponds to 
a known differential equation satisfied by it generating function. In Table [3] we compute the asymptotic 
number of some of these graphs. 

We make one simple observation. The {fc, k + 2}-regular graphs are isomorphic to the k + 2-regular graphs 
with loops, by simply adding loops to the vertices of degree k. This gives a family of identities 



3.2. Set covers and uniform hypergraphs. We now illustrate the method on another family of objects, 
which results in set covers and uniform hypergraphs. An n-set is a set of cardinality n. 

Definition 1 (fc-cover of a set). A collection of r-sets B = {B\, B r } is an r-cover of S if [Jl =1 Bi = S . 
If S = [n] = {1, 2, . . . , n}, then it is an r-cover of n. A cover is restrictive if all of the Bi are distinct. An 
r-cover is fc-regular if any given element occurs in exactly k-subsets. 

A combinatorial argument shows that the number of distinct covers for a set of n elements is 



which is clearly not P-recursive (equivalently, its generating series is not D-finite.) 

Devitt and Jackson [lOj give a generating function for the number of fc-regular r-covers of [n], a notion 
introduced by Comtet [8] . Further, they prove that the number of arithmetic operations required to actually 
calculate the number of fc-covers of an n set by their method is bounded by cn k logn. Results for fixed k, 
specifically k — 2,3 were treated by Comtet [8] and Bender [1] respectively. 

We can derive enumeration formulas. For example, a fc-regular graph on n vertices is a restrictive fc-cover 
of [n] into 2-sets. In general, calculating the generating function for restrictive fc-covers of [n] into j-sets can 
be expressed as 



This yields the following simple consequence of Theorem 12.11 

Corollary 3.2. Let S be a finite set of integers. For fixed n, and fixed k, the exponential generating function 
for k-regular S-covers of sets is D-finite, and is given by the scalar product 



Example. We can express the problem of counting distinct restrictive 2-covers of a set of 
cardinality n by sets of cardinality less than 5 as a scalar product. Denote the exponential 
generating function of such set covers, by S(t). We have, 



(E[e 2 ],exp {t(hi + h l+2 ))) = (E[h 2 ], exp(th i+2 )} . 





To determine fc-covers with mixed-cardinality sets, say both i and j, we calculate 





S(t) = (E[ ei + e 2 + e 3 + e 4 ],exp(th 2 )} . 
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This problem is perfectly suited to either of our algorithms. We can determine this 
differential equation, and the initial terms of the counting sequence: 



1, 0, 1, 8, 80, 1037, 17200, 350682, 8544641, 243758420, 8010360039. 

It is worthwhile to remark that for a fixed j, the set coverings by j-sets are equivalent to loopless j-uniform 
hypergraphs without multiplicities. These are encoded by E[ej]. If we wish to encode hypergraphs with 
loops, we replace by hj, and if we wish to encode hypergraphs with multiplicities we replace E by H. 

4. Asymptotic analysis 

Now that we have established how to determine the differential equations satisfied by regular families of 
combinatorial objects, we process these differential equations to obtain asymptotic enumeration results. 

Asymptotic enumeration of regular graphs is a topic that has received a great deal of attention. Indeed, 
as Gropp [16] points out, the basic problem of regular graph enumeration was considered before graphs were 
even "invented" , over 120 years ago. We first see some explicit results for graphs with fixed degree sequences 
in the work of Read [55] HS]> however, these are rumoured to be "difficult to penetrate". Nonetheless, 
one can determine an asymptotic expression for the number of 3- regular graphs. Bender and Caufield [2] 
produce the first general asymptotic formula for the number of fc-regular graphs on n vertices, and Bollobas 
produces a similar result by a more probabilistic approach that generalizes with ease to treat hypergraphs. 
Next, work by McKay [20] and McKay and Wormald [23] consider the problem of k which is not fixed, 
but rather a function of n, and they achieve a formula which they believe to be true in general, valid as 
n — > oo uniformly for 1 < k = o(n 1 / 2 ) 

(nk)l ( k 2 -l k 3 r2 



(4) 9k(n) ~ wm^W exp I — i — ^ + 0{k ,n) 

This resembles Bollobas' asymptotic formula [4] for labelled fc-regular r-uniform hypergraphs, on n vertices 
when ^ is an integer 

- pt^Pt^ exp - 1)(fc - 1)/2) • 

He also gives a formula for hypergraphs in which hyperedges only have single vertex intersections, which 
gives the constant in the McKay and Wormald formula for r = 2. 

The asymptotic enumeration problem of regular graphs has been treated with a variety of methods, such 
as the multidimensional Cauchy integral technique mentioned earlier [20] , a "switching" technique based on 
inclusion exclusion [211123] . and some direct combinatorial arguments on the equivalent problem of symmetric 
(0,1) matrices with fixed row sum [2]. 

Bollobas [5j remarks that the number of /c-regular unlabelled graphs grows asymptotically like l n /n\ as n 
tends to infinity and where l n is the number of labelled regular graphs. Intuitively, this is due to the fact that, 
for most large graphs with no isolated vertices, and at most one vertex of maximal degree, the automorphism 
group consists of only the identity automorphism. 

The enumeration of other configurations is relatively untreated. Gessel remarked [14] that the exponen- 
tial generating functions of fc-regular r-uniform hypergraphs (with and without loops, with and without 
multiplicities) are D-finite and the differential equations they satisfy are obtainable via the scalar prod- 
uct. Domocos, [11] determines a scalar product form for the generating series of minimal coverings that are 
multipartite hypergraphs. 

Here we continue to treat a variety of configurations. The results are tabulated in Tableland Tabic [4] and 
allow for a comparison across objects rather than regularity parameter. All of the results were automatically 
generated. 

4.1. Technique. Our method is a classical singularity analysis of formal solutions of the linear differential 
equations. It is precisely the same method we used in our analysis of fc-uniform Young tableaux [7J, and 
thus we do not repeat the details here. Instead, after a short description of the major steps, we present the 
fruits of our analyses. A Maple worksheet of the computations is available at 
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k 




1 


2 


3 


Formula 






n\ 

V n 


/3\ f n\{n/2)\ 


I O / (r, 


\2J n 


Simple graphs 
E[e 2 ] 


Sloanc # 
Constant 


A001147 
1 


A001205 


A002829 
.043 


Graphs with loops 
E[h 2 ] 


Sloanc # 
Constant 


A001147 
1 


A108246 

1 . . — 

e 4 /V7r 


A110039 
.318 


Multigraphs 
H[e 2 ] 


Sloanc # 
Constant 


A001147 
1 


A002137 


A108243 
.318 


Multigraphs with 

loops 

H[h 2 ] 


Sloanc # 
Constant 


A001147 
1 


A002135 


A005814 
2.35 



Table 3. Asymptotic enumeration formulas for different classes of /c-regular graphs. For- 
mulas for 1- and 3- regular are valid only for even n. 



In the simplest cases, essentially the cases we could analyze directly with combinatorial arguments, we 
can solve the differential equation and do an asymptotic analysis on the solution. In the more complex 
cases, we first convert our differential equation to the recurrence satisfied by the coefficients. Our series are 
D-finite, and thus such a sequence is bounded by a rational power of n\, and thus, we scale our sequence 
until it is convergent, and this allows us a more precise analysis. We convert this recurrence back to a 
differential equation, and determine the roots of the polynomial which is the coefficient of the leading term. 
From this we can calculate the dominant singularity, and determine a power series solution to the differential 
equation around this point. From this, we analyze the solution to determine an asymptotic expression for the 
coefficients. This can be done automatically using tools from Maple, specifically DEtools and gfun. Finally, 
by generating sufficiently many terms in the sequence, we compare with the formula to determine a value 
for the constant. 

In the tables that follow, the Sloane sequence number refer to the counting sequence as indexed in the 
Sloane On-line Encyclopedia of Integer Sequences [S7J . 

Furthermore, we mention only the dominant term of the asymptotic expression, but we could get the 
subsequent terms, save for the appropriate constants. This is a consequence of the principal weakness of our 
method- we cannot generate exact expressions for the constants. 

We remark that, as presented in Tabled graphs of different types have the same asymptotic development, 
but differ only in the constants. There are a few observations we can make based on this data. We see 
that allowing repetitions influences only the constant of the asymptotic expansion, and often only slightly. 
We see this again in the cycle covers in Table |H The formulas look different than the graph formulas we 
presented earlier, however if we expand the factorials with Stirling's formula for n!, we quickly see that they 
are the same. 

Although we are able to compute the differential equations for the generating functions of the classes of 
4- regular (mufti-) graphs (with loops), their asymptotic analysis is more complicated to do in an automated 
fashion, because a saddle point analysis arises. This is an obvious starting point for future work. Most of 
the tools are already implemented, it is mostly a question of understanding them, and determining how to 
best automate them. 

We could make further observations by considering directed versions of any of these structures. For 
directed graphs, we need only consider E o L2, where is the species of lists of length fc, and we could 
generalize hypergraphs in a number of ways; by putting an order, or even an orientation on each "edge" 
using the species of cycles or lists, as in the directed graph case. 
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fc or S 


Restrictions 


Sloanc # 


Formula 


Constant 


5-regular graphs 


{1,2} 




A00986 


n -?e V2n n\ 


e^/2^ 


{2,3} 




A110040 




0.007 








n- 1 -j n!(n/2)! 




{1,3} 


n = mod 2 


A110039 


0.43 


{1,2,3} 




A110041 




0.05 


fc-regular 3-uniform hyper-graphs: S[es] 


1 


n = mod 3 


A025035 


/1\" /3 n\ 
W («/3)! 


1 


2 


n = mod 3 


AllOlOO 


[jj "K-/3J! 


0.175 


3 




AllOlOl 




0.037 


fc-regular 4-uniform hyper-graphs: E[e^\ 


1 


n = mod 4 


A110102 


/l\" /4 n! 
1,4! J (n/4)! 
/2\ n//2 


1 


2 


n = mod 2 


A110103 


n" 1 (J J n!(n/2)! 


0.100 



Table 4. Asymptotic enumeration of different classes of regular objects 



5. Comments, conclusions and perspectives 

5.1. Asymptotic expansions of different families of functions. Coefficients of taylor expansions of 
algebraic functions have a known kind of expansion, that can in fact we used to establish the transcendence 
of some series [12] . We can also describe asymptotic discrepancy criteria for coefficients of D-finite functions. 

As we remarked earlier, coefficients of D-finite series are also restricted in their asymptotic growth. A 
more complete version of this criteria is following theorem presented in Wimp and Zeilberger [30j . 

Theorem 5.1 (Wimp and Zeilberger). Suppose that f(t) = ^ n>0 fnf 1 is a D-finite series in in C[[t}]. 
Then, for sufficiently large n, the coefficients f n have an asymptotic expansion which is a sum of terms of 
the form 

A(n!) r / s exp(Q(n 1 /" i ))w"n Q (logn) fc , 
where r, s, to, k G N, Q is a polynomial and A, oj, a, are complex numbers. 

We may ask ourselves, have our examples encompassed the full asymptotic potential of D-finite sequences? 
We are very curious about the combinatorial structure of families which do have such expansions. Are D-finite 
species sufficient to consider? 

For example, in our earlier study of fc- uniform Young Tableaux [7], which are enumerated by the scalar 
product (H[e\ + e2\ 1 exp(hkt)) have a conjectured form (verified for fc = 1..4) of 

Vn V2 I 2tt ) U - { fc! J n fe / 4 ' U 

This is an exact conjecture, more complete than the examples we have presented here, although it results from 
the same kind of calculation, and presumably if we completed the complex saddle point analyses required 
for the fc = 4 cases, we might be able to guess such a form for our examples. 

It is also of interest to note that while in some cases operations of summation and integration preserve 
D-finiteness, G r , the class of all regular graphs is not D-finite. The same is true of all the classes we have 
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presented here: Although for any k, the subclass of /c-regular objects is D-finite, the larger subclass of regular 
objects is not. This is interesting and can help us refine our notion of D-finitcncss. 

5.2. D-finite species? We have thus far restrained ourselves from defining a notion of D-finite species. 
Ideally, such a theory would contain two main components: A "D-finite species" should satisfy some sort of 
system combinatorial differential equations with polynomial coefficients; and symmetric series of such species 
should be D-finite. We would then expect to be able to have theorems of the form: 

(1) If F and G are D-finite species, then so are F + G, and F • G, F'; 

(2) If G is a polynomial species, (in particular, if its cycle index series is a polynomial), then F o G is a 
D-finite species; 

(3) If F satisfies an "algebraic equation" of species, including for instance, equations of the form F = 
XP(X, F) for polynomial species P, then F is a D-finite species. 

A candidate definition is given in [24] . however more work remains to be done. 

Finally, much work has been done to characterize combinatorial classes of objects with rational and 
algebraic generating series (see [5] for a recent summary), and hopefully this work is a step towards such 
a characterization for D-finite generating functions. We are encouraged by a recent thought of Flajolet, 
Gerhold and Salvy [T5] . 

Almost every thing is non-holonomic unless it is holonomic by design. 

(Series are holonomic if and only if they are D-finite.) They follow this with the remark that there are several 
surprising exceptions to this rule, notably fc-regular graphs. Hopefully we have demonstrated that this isn't 
so surprising; That in fact, there are deep reasons underlying the D-finiteness of objects with this sort of 
regularity, and that furthermore, this D-finiteness can be exploited in an automatic way. 

Acknowledgements. The author wishes to thank Francois Bergeron, Bruno Salvy, Frederic Chyzak, and 
indirectly Frederic Jouhet, for fruitful, instructive discussions. Thank you also to Cedric Chauve for useful 
comments on the text. 

Sequences: A001147, A001205, A002829, A108246, A110039, A002137, A108243, A002135, A005814,A00986 
A110040, A110041, A025035 A110100,A110101, A110102, A110103. 
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